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The notion of fidelity in quantum information science has been recently applied to analyze quan- 
tum phase transitions from the viewpoint of the ground-state (GS) overlap for various many-body 
systems. In this work, we unveil the intrinsic relation between the GS fidelity and the derivatives 
of GS energy and find that they play equivalent role in identifying the quantum phase transition. 
The general connection between the two approaches enables us to understand the different singu- 
larity and scaling behaviors of fidelity exhibited in various systems on general grounds. Our general 
conclusions are illustrated via several quantum spin models which exhibit different kinds of QPTs. 
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I. INTRODUCTION 



Quantum phase transitions (QPTs) that happened at 
the zero temperature is purely a phenomenon of ground 
state (GS) transition driven by external parameters. Tra- 
ditionally, QPTs are described in terms of order parame- 
ter and symmetry breaking within the Landau-Ginzburg 
paradigm which have been extensively studied in con- 
densed matter physics Q. In recent years, QPT has 
also attracted a lot of attention in quantum-information 
science 0, in which one of the research focus is the 
role of quantum entanglement in characterizing QPTs 
[3, H, 0, H, @] ■ More recently, another concept in quan- 
tum information science, i.e., the fidelity has been put 
forward to identify QPTs from the perspective of the GS 
wave functions [3, 0] ■ The GS fidelity is defined as the 
overlap between two ground states with only slightly dif- 
ferent values of the external parameters [8] and thus is a 
pure geometrical quantity. Since no a priori knowledge 
of the order parameter is needed, the fidelity might be a 
potential universal criteria for characterizing the QPTs 
E i, [M E, [il Q m • An increasing interest has 
been drawn in the role of GS fidelity in detecting QPTs 
for various many-body systems d, [ifl, [HI, [13, [H, 113, [lB| , 
since Zanardi and Paunkovic first exploited it to iden- 
tify QPTs in the XY spin chain [s*] where the fidelity 
shows a narrow drop at the transition point. Remark- 
ably, the success of fidelity analysis in dealing with the 
Bose-Hubbard model [5| and spin systems [3, [lB| implies 
that it may have practical relevance even for more com- 
plicate strongly interacting systems where no a simple 
description is possible. The relation between the fidelity 
and Berry phase [l^ has also been unveiled in terms of 
Riemannian metric tensors • 

Generally one may expect that the GS fidelity shows 
an abrupt drop in the vicinity of the QPT point of the 
system as a consequence of the dramatic change of the 
structure of the GS. This is true for the first-order QPTs 
caused by a level crossing of GSs due to the GS wave- 
functions at the different sides of the level-crossing point 
are almost orthogonal (orthogonal at the level-crossing 



point). However, this is neither an obvious nor a gen- 
eral conclusion for a continuous QPT where no GS level 
crossing occurs and the GS evolves "adiabatically" in 
the parameter space [lll[i3|- Conventionally, QPTs are 
characterized by singularities of the ground state energy: 
first-order QPTs are characterized by discontinuities in 
the first derivative of the energy, whereas second-order 
(nth-order) QPTs are characterized by discontinuities in 
the second (nth) derivative of the energy. An impor- 
tant question is whether the singularity of GS energy is 
intrinsically related to the GS overlap? Answering this 
question is no doubt significant for a deeper understand- 
ing of QPTs, and the validity and limitation of fidelity 
as a measure of QPTs. 

So far the studies of fidelity as a measure of QPTs are 
based on the analysis of particular many-body models 
[1, [§] . The general connection between GS fidelity and 
QPTs is not yet well established. In this work, we shall 
discuss, in a general framework, how the GS fidelity can 
be related to a QPT characterized by non-analyticities of 
the GS energy. Our result shows that the singularity and 
scaling behaviors of the GS fidelity (or it derivatives) are 
directly related to its correspondences of the derivative 
of GS energy. 

Our paper is organized as follows. In section II, 
we reveal the general relation between the fidelity sus- 
ceptibility and the 2nd order derivative of the GS en- 
ergy. The subsequent sections are devoted to two exam- 
ples which exhibit the 2nd order and Kosterlitz-Thouless 
(KT) QPTs, respectively. A summary is given in the last 
section. 



II. FIDELITY AND QUANTUM PHASE 
TRANSITION 

The general Hamiltonian of a quantum many-body sys- 
tem undergoing QPTs reads 



H{X) = Ho + \Hi, 



(1) 



where Hi is supposed to be the driving term with A the 
control parameter. In terms of the eigenstates |5'„(A)) of 



2 



H{X), the Hamiltonian can be reformulated as H{X) = 
-Bn(A)l*„(A))(^'„(A)|, where Af is dimensions of 
the Hilbert space. The GS fidelity is defined as the over- 
lap between |^'o(A)) and |*o(A + 6)), i.e. 



F{X,6) = |(vI/o(A)|*o(A + (5)}|, 



(2) 



where (A) is the GS wavefunction corresponding to the 
parameter A and (5 is a small quantity Q. It is obvious 
that the fidelity is dependent of 5. The rate of change of 
fidelity is given by the second derivative of fidelity^] or 
fidelity susceptibility (FS) [IH 



5(A) = a|F(A,5) \s^o- 
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(52 



(3) 



which is S independent and sometimes a more effective 
quantity to detect the QPT. 

First we discuss the first order QPT (IQPT) which is 
induced by GS level crossing. We assume Eo{X) < Ei{X) 
for A < Ac with Ac being the crossing point where the 
GS energy level Eq crosses over the first excited level 
El, therefore we must have Ei{X) < Eo{X) for A > Ac- 
Explicitly, for the IQPT, the GS energy is defined as 



Eg{X) 



EoiX) A < Ac 
^i(A) A>Ac 



The derivative of the GS energy is defined as 



dEgjX) 
dX 



lim 



£o(A)-go(A,) 

x^K-a A^ri: 

ImiA^A.+O A-A. 



_ dEqjX) . . 

_ ^E^{\) . . 
- 9A > Ac 



In general, -^Eq (A) is not equal to -^Ei (A) , there- 
fore the first derivative of GS energy for a IQPT 
is not continuous at the transition point. Accord- 
ingly, the GS fidelity at the transition point is de- 
fined as F{Xc,S) = |(*g(Ac- V2)|*3(Ac + <5/2))| - 
(*o(Ac - (5/2)|^'i(Ac (5/2))| , therefore we have 
lim^^o ^(Ac, <5) — which means that there appears a 
sharp drop in the transition point. It is straightforward 
that the sudden drop in the transition point shares the 
same physical origin with the discontinuity of the first 
derivative of GS energy. 

Next we consider the continuous phase transition for 
which the GS of the Hamiltonian is nondegenerate for a 
finite system. For the case of 2nd order QPT, there is no 
level crossing for the GS energy, therefore we can always 
represent Eg (A) = Eq (A). The first derivative of the GS 
energy is given by ^Eq (A) = («'o(A)| Hi |«'o(A)) , which 
is noting else but the Hellmann-Feynman theorem. It 
is straightforward to get the second derivative of GS 
energy ^EoiX) = («'o(A)| ffi |9a*o(A)) + h.c. . In- 
serting the identity operator J2n=o l*n{A)) (^n(A)| = 1 
between Hi and |9a*o(A)}, we then get ■^Eo{X) — 

En=iHL{'^niX)\dxM^)) + h.c. With i7i„ = 
(*o(A)| Hi |*„(A)) . We note that the term of n = is 
not included in the summation because it gives zero due 



to (9a*o(A) |*o(A)) + (*o(A) |9a*o(A)) = 0. Differen- 
tiating the eigenequation 7?|^o(A)) — i?o |^o(A)) with 
respect to A and taking the inner product with 5'„(A) 
yields (*„(A)|c'aH|*o(A)) + (4'„(A)| H |aA*o(A)) = 
dxEoSno + Eo i^ifniX) |9a*o(A)) . Exploiting the 
hermiticity of H to write (*„(A)| iJ |5a*o(A)) = 
(l'n(A) |<9a*o(A)) , it follows that for n ^ 0, 
(^'„(A) |aA*o(A)) = (*„(A)|9Ai/|*o(A))/(^;o-^^»). 
Therefore, we have 



(9A2 



71^0 



(^„(A)|gi|vI'o(A))|^ 
So(A)-^„(A) 



(4) 



To see the relation of FS with the second derivative of 
GS energy, we rederive the expression of FS [HI, [13] by 
expanding the wavefunction |4'o(A -I- (5)) in the basis of 
eigenstates corresponding to the parameter A, to the first 
order, which leads to 



|*o(A + <5)} 



|^o(A))+(5 



^ gio(A)|vl/„(A)) 
2-.^ EoiX) - E^iX) 



where c = ^1 + 6^ E„^o l^no (A)l' /[EoiX) - E^XW is 
the normalization constant. It follows directly FiX,6) = 
c. Substituting it back into Eq. we get 

5(A) = > '^^^^^1^^^^^^^^^::^. (5) 



|( vl/„(A)|gi|M>o(A) )r 
[EoiX) - E„ iX)Y 



Comparing Eqs. ^ and ([U, we find that for both of 
them the singularities come from the vanishing energy 
gap in the thermodynamic limit. However, a divergence 
does not necessarily occur in the critical point, for ex- 
ample, for the KT QPT as we shall discuss later, no 
singularity arises even the energy gap tends to zero in 
the thermodynamic limit as matrix elements H^q may 
also vanish simultaneously. Despite the apparent simi- 
larity between Eqs. ([5]) and ([3]), we note that fidelity 
susceptibility cannot be expressed by the GS energy 2nd 
derivatives alone [l^. Nevertheless, the different power 
in the denominators in Eqs. ([5]) and (|4]) shows that fi- 
delity susceptibility might be a more sensitive tool to de- 
tect critical points because it can be more singular than 
the second derivative of GS energy. Therefore, for the 
case where the GS energy 2nd derivative is divergent at 
the critical point, it is no doubt that the fidelity suscep- 
tibility is divergent too. However, for the case where the 
second derivative of GS energy is not divergent, (for ex- 
ample the 3rd-order QPTs and the KT phase transition) , 
it is hard to judge whether the fidelity susceptibility be- 
comes singular or not at these critical points just from 
the relation between Eqs. ([5]) and (|4|). 



III. TRANSVERSE FIELD ISING MODEL 

To understand the equivalence of the two approaches 
in characterizing the quantum criticality, we first apply 
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them to study a specific example, say, the transverse field 
Ising model with the Hamiltonian given by 



HiX) 



i=l 



(6) 



where A represents external magnetic field along the z 
axis and the periodic boundary is assumed. Following the 
standard procedure the model can be diagonalized 

by applying Jordan- Wigner transformation which maps 
the spin model to a fermion model and then using the 



Bogoliubov transformation = cos ^ 
The result is given hy H — J^k^k (cl 



Ck 



1 1 with LJk = 



,2 27rfe 



— I (a — cos , where the summation is over 
k = —M, ■ ■ ■ , AI with the assumption of i = 2M + 1. 
The ground state of H is the vacuum \g) of Ck given 

tM 



by \9) 



n 



fe=l 



COS 



with the ground state energy Eq (A) = ~J2k^k, where 
|0)^. and are, respectively, the vacuum and sin- 

gle excitation of the kth mode bk and 9k is defined 
by sin 9k = sin ^jj^ /ujk- It follows that the fidelity is 



nven by F(A, 5) 



11/0=1 2 



cos ■ 



with 9k = 9kiX + S) 



[8|. After some very simple algebras, one get <S'(A) = 
2i-F(A^ = J YXLi {^f ■ In the thermodynamic limit 
with L oo, the summation can be replaced by the in- 



tegral and thus we get 



eo(A) 



1 

2^ 



■ 2 

sm Lp 



^ \siv? ip + {\ — cos if)'' 



i3/2 



and 



167r 



sin^ (/? 



-TT [sin Lp+{\ — cosLpY 



J dip 



where eo (A) = Eq (A) /L and s(A) = S{X)/L. Both of 
the above integrals are divergent when A = 1. Their 
singular behavior can be analyzed in the vicinity of criti- 
cal point Ac = 1 with the asymptotic behavior described 
by -d^eo (A) /d\^ ^ -0.31187 In | A - Ad + const and 
Ins(A) = — In |A — Ad + const, respectively. 

To get further insight on how the two approaches work, 
we study both the scaling behaviors of FS and second 
derivative of GS energy. As shown in Fig. 1, we display 
the second derivatives of GS energy for different lattice 
sizes. Despite no real divergence for finite L, it is obvi- 
ous that the curves exhibit marked anomalies with height 
of peak increasing with the lattice size. The parameter 
Am labelling the position of peak approaches the critical 
point Ac = 1 in a way of Am = 1 — constL~^'^'^'^ . The 
value of peak diverges logarithmically with increasing lat- 
tice size as — 9^eo(A, L)/9A^|a,„ — 0.31132 InL -I- consi. 
By proper scaling scheme, we can fit all the data of 
F = [l - exp((9^eo/i9A^|A„ - 9^eo/9A^)] as a function of 
L (A — Am) for different L into a single curve as displayed 
in the inset of Fig. 1. Correspondingly, the FS exhibits 




FIG. 1: (color online) The derivatives d^Eo (A) /d\^ as a func- 
tion of A. 
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FIG. 2: (color online) The fidelity susceptibility as a function 
of A. 



similar scaling behaviors as shown in Fig. 2. Around the 
critical point, all the data for different L collapse into 
a single curve of F = [1 — s(A)/s(Am)] as a function of 
i (A — Am)- A related quantity has been used in [l^ to 
make scaling analysis for the quantum Ising model. It is 
clear that d'^eo (A) /dX^ and s(A) exhibit similar critical 
behavior around the critical point. 



IV. XXZ SPIN CHAIN 

Having demonstrated the equivalent role of the two dif- 
ferent approaches in identifying QPT on the transverse 
Ising model, we now use this procedure to understand 
the more elusive KT-type QPT for which, contrary to the 
previous expectation, the GS fidelity and the FS turn out 
to be insensitive to QPT in some cases, for example, in 
the spin Ji — J2 model [l^. We now consider the XXZ 
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FIG. 3: (color online) The fidelity susceptibility and second 
order derivative of GS energy for XX Z model as a function 
of A. The data labelled by x and + are obtained by ED and 
they comply with the data by Bethe-ansatz method exactly. 



model which provides a further test for the intrinsic re- 
lationship between the derivatives of GS energy and FS. 
The model is defined by 



H{\) 



1 ^ 



i+1 



i+l) 



(7) 



with A the exchange anisotropy parameter and it is re- 
lated to —H{—X) by a unitary transformation. It is 
well known that the model is in a critical phase for 
— 1 < A < 1, an antiferromagnetic phase at A < —1, 
and ferromagnetic phase at A > 1. The QPT at the fer- 
romagnetic isotropic point A = 1 is a IQPT caused by 
the level crossing of GS. Here we shall focus on the QPT 
at the antiferromagnetic isotropic point A = — 1 which is 
known to be of KT type. In the whole parameter regime 
the model is integrable. The eigenstate of H{\) has the 

a{ni, ■■■ , nM)S-^ ■ ■ ■ S'lF) 



form of 1 4') = V 

with the coefficients having the Bethe-ansatz form [20| 



ni<---<njv/ 



a{ni, ■ 



a</3 



where P is any permutation of (ni, 



phase 0ip,q) = 2tan ^ ( eos[(p+g)/2]~Acos[(p-g)/2] ^ 

eigenenergies E = —-127-1-2 ^-(A — cos ki) are determined 
by a set of quasi- momentum pi {i = 1, • ■ • , M) for M 
down spins which arc the solution of Bcthc ansatz equa- 
tions (BAEs) Nki = 2iTli — X]j=i ^(^i; where li are 
integer or half-odd integer. 

By numerically solving the BAEs, we can get the GS 
state energy for different sizes of L. As shown in Fig. 3, we 
calculate the 2nd derivative of GS energy by both Bethe- 
ansatz and exact diagonaliztion (ED) methods and the 
FS by ED. Our numerical results show that there is no 
any singularity around the critical point Ac = — 1 for 
both the FS and 2nd order derivative of GS energy for 



Asin[(p-g)/2] 



and the 



The 



different sizes. For the GS energy derivatives, this is 
true even in the thermodynamic limit where actually, as 
proven in [21], all the nth order derivatives of GS en- 
ergy are continuous. Thus, in view of its connection to 
the derivatives of GS energy, we can understand why the 
GS fidelity failed to reproduce quantum critical behav- 
ior by the desired singularities in the critical point for 
the KT transition. It is also natural to explain why 
the finite size scaling of FS has a clear distinction be- 
tween noncritical and critical phases because the GS en- 
ergies are known to fulfil quite different finite size scal- 
ing in critical or noncritical phase [13, [2^ . The GS en- 
ergy in the critical regime fulfills the finite-size correction 
£:o(A, L)/L = eo(A) - ttu (A) c/6L^ + o{L-^) with c = 1 
being the conformal anomaly number and v — tt sin 7/7 
the "sound velocity" for the model 22]. Here 7 is defined 
by A = — cos 7. On the other hand, the conformal invari- 
ance is broken for A < —1 and the GS energy fulfills a 
different finite size scaling [23] . Bearing in mind that the 
equivalent role of FS and derivative of GS energy, we ex- 
pect that the FS exhibits similar scaling behavior in the 
critical regime as its counterpart of derivative of GS en- 
ergy does and there is a deviation from the scaling behav- 
ior when A < — 1. Indeed, the FS has been shown to scale 
differently in critical and noncritical phases based on the 
results by ED method T?]. This feature can be used as 
an indicator of KT QPT occurring in this model. Before 
ending the discussion, we would like to give a remark to 
the problem whether the fidelity approach is suitable to 
detect the KT transition? As no rigorous result of the 
FS even for the exactly solvable model is available in the 
thermodynamic limit, therefore this question is still in a 
doubtful status [3, [l§| and remains an open problem. 



V. CONCLUSIONS 

We have shown the intrinsic relationship of the fidelity 
and derivatives of GS energy and revealed their equiva- 
lent role in identifying the QPTs. Within our framework, 
the singularities and scaling behavior of fidelity for vari- 
ous systems can be well understood on general grounds. 
Through concrete examples, we display that divergence 
of the FS or derivatives of GS energy can be applied to 
identify the 2nd order QPT, whereas for the KT transi- 
tion the criticality is not a sufficient condition to ensure 
divergence of the FS or derivatives of GS energy. 
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